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The lineshape of the tunnehng conductance in double quantum wells with a large-scale rough- 
ness of heterointerfaces is investigated. Large-scale variations of coupled energy levels and scattering 
VO \ due to the short-range potential are taken into account. The interplay between the inhomogeneous 

. broadening, induced by the non-screened part of large-scale potential, and the homogeneous broad- 

ening due to the scattering by short-range potentials is considered. It is shown that the large 
inhomogeneous broadening can be strongly modified by nonlocal effects involved in the proposed 
■ mechanism of inhomogeneity. Related change of lineshape of the resonant tunneling conductance 

(~| ' between Gaussian and Lorentzian peaks is described. The theoretical results agree quite well with 

, experimental data. 

CD 



I. INTRODUCTION 



. Resonant tunneling in semiconductor heterostructures has been widely investigated ever since Tsu and Esaki pro- 
Q ' posed the double-barrier resonant-tunneling diode (see Ref. Q for a recent review). New developments came 
O , through from studies of interlayer tunneling spectroscopy between parallel two-dimensional electron systems (2DES) 
using the technique of independent contacts to closely located 2DES. ||,^. The 2DES arc formed in two GaAs quan- 
^-H \ turn wells (QW) separated by a Al^^Gai-^^As barrier. Because the in-plane momentum and the energy are conserved, 
^ • the 2D-2D tunneling current exhibits sharp resonance peak whose broadening is determined by different collision pro- 
■ cesses in the nonideal double quantum well (DQW) structure. This property allows to study scattering mechanisms 
through tunneling spectroscopy method. |^ Furthermore, broadening effects may be important in a novel quantum 
transistor based on 2D-2D tunneling in independently contacted DQWs. ^ 

The aim of this paper is to describe the lineshape of the resonant tunneling current in nonideal DQWs with 
• independent contacts to each QW, when, in addition to usual homogeneous broadening induced by short-range 
, scattering, the inhomogeneous broadening due to large-scale variations of heterointerfaces is taken into account. The 
^ ■ latter scattering mechanism has an essential effect on the form of the peak, because smooth variations of the DQW 
Ci ' energy levels due to large-scale random variations of the widths of right (r-) and left QWs can not be screened, 
^ , even though the screening potential involves all possible redistributions of electrons within the DQW structure. 
I ■ Even though the averaged large-scale potential is screened in heavily doped structures, the intersubband energy is 
still nonuniform over the plane of the quantum well. In Fig. 1, a schematic view of the band diagram of DQWs and 
spatial variations of the energy levels are depicted for illustration. Our theory is valid when the DQW width is smaller 
^ ' than the correlation length £c for nonuniformities of the heterointerfaces in the DQW. A very similar mechanism was 
recently proposed in a single quantum well for describing the inhomogeneous broadening of intersubband transitions, 
. ^ with one subband occupancy, Q and for new effects in classical magnetotransport in the case of double subband 
^ . occupancyJI 

?H We show that the Lorentzian lineshape for the tunneling current peak, in the case of short-range collision-induced 

broadening, assumes a Gaussian shape due to the inhomogeneous broadening, if nonlocal effects are discarded due 
to sufRciently large ic- However, for not too large £c, we obtain the transformation from a Gaussian to Lorentzian 
lineshape due to nonlocal effects on the inhomogeneous broadening. Moreover, inhomogeneous and nonlocal effects 
essentially modify the half width at half maximum (HWHM) of the peak. As it is shown below, our theoretical results 
are in quite reasonable agreement with the experimental ones of Ref. . 

The paper is organized in a following way. In Sec. II we evaluate the expression for the tunneling current up 
to second order in the weak interwell tunneling coupling and use the path-integral representation to calculate the 
tunneling conductance in terms of the averaged product of Green's functions for electron in left (l-) and right (r-) 
QWs. The lineshape of the resonant tunneling conductance is analyzed in Sec. Ill in a quasiclassical approximation. 
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The list of assumptions and concluding remarks are given in the Sec. IV. The Appendix A contains estimates of 
the parameters used in the nonscreened potential, due to large scale nonuniformities of the heterointerfaces, and in 
Appendix B we briefly discuss the optimal fluctuation method and the straightforward trajectory approximation used 
in Sec. III. 



II. TUNNELING CURRENT 



Electron states in I- and r-QWs are described by the Hamiltonians 

Zm 



P 

2m 



Ur 



(1) 



where A is the interlevel splitting without tunneling and m is the effective mass. The effect of fluctuations of 
heterointerfaces and scattering processes are described by large-scale and short-range potentials Ui^rx and Ui^rx in I, 
r-QWs. The screening potential, Vx, included in -ff;,r, is determined from the Poisson equation (see Appendix A) and 
only the averaged large-scale potential is screened as 



U 



Ur 



+ Vx = 0. 



(2) 



Taking into account the interwcU tunneling coupling, we use a 2 x 2 one-electron Hamiltonian matrix as 

hi T 

T hr 



(3) 



where the diagonal terms are given by 



hi 

hr 



2m 



SUx 
2 



Urx r, 

2m 2 



the non-screening part of the large-scale potential is SUx = Ui^ — Urx and the nondiagonal terms are given by the 
tunneling matrix element T (the coupling energy) . In the following, we assume that the random potentials introduced 
above are statistically independent and described by Gaussian correlation functions 



(C/,xC//x) = S,/W,i\^ - x^, {U,xU,^ = S,/Wji\^ - x'l), 

{SUxSUx') = Wli\^ - x'l) + Wr{\^ - x'l) EE W(|X - X' 



(4) 



where the functions Wi^r{^) and w{x.) are discussed in Appendix A. We also neglect here the in-plane variations of 
the matrix element T (see discussion in Ref. |^). 

The interwell tunneling current is expressed in terms of the density matrix pt according to |P,|l0| 



J± 



\e\T2_ 

h £2 



tT{(7yPt), where pt 



Pit Pt 

pt Prt 



(5) 



with ay being the y component of the Pauli matrix and the trace includes both the average over large-scale and 
short-range random potentials and the summation over electron states. Non-diagonal and diagonal components of 
the density matrix in Eq. (||) are connected by the relation {6 +0) 



Pt = '^ I dt'e'''e-^^^('-''y\pu>-prt')e 



ihr(t-t')/h 



(6) 



Using a set of wave functions (x|j'A) = tA^x which are determined by the eigenvalue problems in the jth QW hjip^^ 



Sxipj^ we rewrite the tunneling current (W as 



2 



Here ((. . .)) means the average over short-range and large-scale potentials. After substitution of Eq. (||) in Eq. (|^), 
we obtain 



L2 



where the above second equation is written for the zero-temperature case and e^j is the quasi-Fermi level in the jth 
QW. 

In order to calculate J±, it is convenient to use the retarded (R) Green's functions for the electron in /- and r-QWs, 
which are defined as 

A ^^^^ -e-iS) 



and the advanced (A) Green's functions given by C7^g(x, x') = tj"g(x',x)*. The tunneling current assumes the form 

r^de jd^l dx' ^ {-if{{giMX)QU^.^))). (10) 



2'Kfl L2 



where k = 1 for a — b and fc = for a ^ b. For small applied voltages satisfying \epi — e^rl ^ Syr,; — e^, we introduce 
the tunneling conductance, G(A), through the relation J± = G{A)V. The interwell voltage, V, is connected with 
the quasi-Fermi level difference by the relation V = {Spi — epr)/s- Then from Eq. ( p^ it follows that the tunneling 
conductance can be written as 

G(A) = 1^x1 dx' E (-l)'=(fe(x,x')e.(x',x))). (11) 

ab—RA 

Furthermore according to Eq. (^, the short-range potentials in the l-QW and r-QW are statistically independent, 
then the two-particle correlation function ((...)) in Eq. ( pi] ) can be rewritten exactly in terms of the Green's functions 
G°g(x, x') — {Qj^{x,x')) averaged over the short-range potentials. The Dyson equation for this Green's functions is 
written as 

{hj - e)G^%(x,x') + J dxiI]^\(x,xi)G^\(xi,x') = ^(x- x'). (12) 

Here the Hamiltonians hi^r coincide with those given in Eq. (??), without the short-range potentials Ui^rx^ and 
S°g(x, x') is the self-energy function. For (5-correlated potentials, we have to use I]"^(x, x') oc (5(x — x'). Neglecting 
the renormalization of energy spectra, we rewrite Eq. (|l2|) in terms of the broadening energy 7j of the jth QW as 

{hj - e T *7,)G;^"(x,x') = 5(x - x'), (13) 

where the upper sign corresponds to and the lower one to G^. 

It is convenient to write the Green's functions through path integrals as |ll[ 



Gf,(x,x') = J / die-'(^+^^'-^)*/'' r %{x,}exp / drimi^l - 5U^^) 

G«(x,x') = J r die-'^^+^T-)*/'' / ' I?{x,} exp [-^ / ^^("ix^ + ,5C/,jl , (14) 



and G^g(x, x') = G^g(x', x)*. The average over the non-screened large-scale potential in Eq. (|l l|) , for a Gaussian- type 
random potential 5C/x, is performed using the well-known exact formula 
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exp / dxf^SU^ ) = cxp 



i y dx y dx7,u;(|x-x'|)/. 



(15) 



for some arbitrary function /x. Since random potentials are involved in both path integrals, we choose these functions 
as 

f^ = ±{i/2h) [ dri(5(x-x^J ± (i/2?i) / dT2(5(x-x^J. 



Using these transformations in the correlation functions of Eq. (njl) , we finally obtain 



^ (-l)^-(G?,^(x,x')G^,,(x',x)) [ diie(^'+^^"^/'' [ dhe''^'^/' 

I, „ . fl J —OQ J — 



lb—RA 



-x' 

_j_g— ieF(ii — t2)//i 



yt2=x 



V{yr} exp [-5+(tit2|Xr, Vr)] 



yo=x 



^'{yrlcxp [-5_(ti<2|Xr,y^)] } +C.C 



(16) 



where the two-particle actions 5±(tit2|xr,yr) are written in the form 



5±(tit2|Xr,yr) 



im 







1 /■ ^ /■ ^ 

Oft Jo Jo 

+ ^ / / cfr'^dy^ -y^-l) 

o/l Jo Jo 

T772 / d,T dr'^dx^ - y^-l). 
4?i Jo Jo 



(17) 



Substituting Eq. ( |16D into Eq. ( 11 ) and making convenient change of variables (in particular, separating the straight 
path according to x,- — > [ur/ti +Xr] and jr [u(t2 — T)/t2 +yT], for integral from exp(— 5+), or y^- [ur/i2 +yT], 
for integral from exp(— we can express G(A) in terms of contour integrals as 



G(A) 



V{xr} f T>{yr}J2\ 6"'^^'*^=^*^^/'* exp 



2h 



X exp 



'2h 



dT±'^± / dryi 



exp [-K±{ti,t2,Xr,yr)] 2 ^ + c.c 



(18) 



where u = x — x'. The contributions of non-screened potentials to the correlation function is given by the factors 

if±(fl,f2,Xr,yr) = — 2 / d'T dT'w{\Xr -Xr' +u{t -T')/ti\) 

oh Jo Jo 

+ ^ I' dr r dT'w{\yr - yr' ± u(t - t')A2|) 
on Jo Jo 

T772 / dT f dT'w±{\Xr-yr' +u{T/ti±T'/t2)\) . (19) 

47i Jo Jo 



with ti;_(|z|) — w{\z\) and w+dzj) ~ w{\z — u|). Note that /v+ comes from averaging both retarded or both advanced 
Green's functions while corresponds to averaging the product of retarded and advanced Green's functions. 
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III. LINESHAPE OF THE CONDUCTANCE PEAK 



In order to calculate the path integrals in Eq. ([l^), we wiU neglect in Eqs. (|l|) deviations Xt and yr in 

the arguments of the correlation function w(|...|) by supposing that these deviations are smaller than £c, i.e., using 
the approach of straightforward trajectory in Eq. (|l9|). We justify such an approximation in Appendix B, where 
the optimal fluctuation method [ p^ is used, in order to extract the optimal trajectories which give the maximal 
contribution to the path integrals. With this approximation, we can calculate the path integrals for the free motion 
exactly and the conductance, given by Eq. (Eq), is rewritten as 



G(A) 



xE 



du 



-oo J — oo ^1^2 



2h 



(20) 



where the factors K± are reduced to 

K±{ti,t2,u) 



1 

8^ 
1 



dr 



dr 



dT'w{u\T - T'\/tl) 

dT'w{u\T — T'\/t2) 



T772 /'^^ r dT'w±{\n{T/tl±T' /t2)\) 

-ih-" Jo Jo 



(21) 



Let us for a moment ignore the contribution from the terms with the upper sign in Eq. (po|). Defining new variables 
X — T lt\^2 and x' — t' /ti^2 in the factor K^iti, t2, u), we obtain the conductance in the form 



G(A) = 



du 



dt2e 



2h 



ft7^ - t,-M 



tlt2 
itl+t2f 



8h^ 



-w 



+ c.c. , 



(22) 



where the large-scale correlation function is transformed as W{u/£c) = Se dx dx' e:cp[—{u/lc)'^(x ~ x')'^] and can 
be rewritten as 

W{x)0 = V^x-i erf(x) - x~^[l - e""'], (23) 
and erf (a;) is the error function. Introducing new time variables t = ti — t2 and t = {ti + 12) jl it follows that 



G(A) 



du 



dt 



^g.[A(t+r/2)-e^r]Ag^p 



i2 -r2/4 

■W 



2h t^-T^/A 2h- 



+ c.c. , 



(24) 



where 7 = 7; + 7r and A7 = (7; — 7r)/2 are the total collision-induced broadening and the broadening difference in 
I- and r-QWs respectively. Since the time scale of r is of the order of h/ep and a typical t is of the order of ?i/7e// 
in the integrals of Eq. (|2J), we can replace t"^ — t"^ by t^, due to the quasiclassical condition ^^ff <C Sp and the 
integration over r gives us 2'!rhS[ep — m{u/t)'^ /2]. After straightforward integration over u we finally obtain 



G(A)=,(f 



exp 



2r 



Vpt 



c.c. 



(25) 



where = m/Trh? is the 2D density of states, the correlation function is given by Eq. (p3|), and Vp is the Fermi 
velocity. 

Consider first the limiting case of the local response, assuming 



(w^n/7,„4)' < 1. 



(26) 
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where the effective HWHM due to both contribution from collision processes and inhomogeneous broadening is 
determined by G{j^ff) = G(0)/2. Under such a condition the correlation function (^) assumes the form W{u/£c) ~ 

W{0) = Se and for conductance lineshape from Eq. (EST) it follows that flj] 



GiA) = 2(^\ /° dte^*/''-(*^*/^'*)'cos[(A/?i)t] 



/ — OC 



(eTT r27/(A2 + 72), _ fc«7 

h \ (V2^/&)exp[-(A/^/2fc)2] &»7 



(27) 



where the limiting cases determine the Lorentzian or Gaussian lineshape. Notice that the Lorentzian shape of the 
peak tails is always found for big enough |A|. 

Now, consider the variables ti,2 in the integral with the upper sign in Eq. (20). They are estimated of the order of 
Ti/sp. Thus u is of the order of Ti/ ^mep <^ £c, and as a result, we can replace the factor exp[— _ft'_|_(ti, t2, u)] by the 
expression exp[— (ti — t2)^w{0) / {Sfi^)] where the exponential factor is of the order of {Se/sf)^ <C 1 and the expression 
can be approximated by the unity. After straightforward integrations over u and r, we are left with an integral over 
t given by 

(eT)2p2D f° ,_^j^,+y^+iA)t/h^-i2ep.t/n , {eT)^P2D 



I- 



Such a contribution can be discarded in comparison with the results from Eqs. ( p5[ ) and (^7|) because this term leads 
to corrections of the order of ^^jj/ep. 

In Fig. 2, we plot G(A), using Eq. (p7[), as function of A/7 for different relative contributions of the scattering 
processes and of the inhomogeneous broadening, i.e., for different ratios fc/7. The solid, dashed, dotted, dot-dashed, 
and dot-dot-dashed curves in Fig. 2 correspond to fc/7 — 0.3, 0.6, 1, 3, and 6, respectively. The factor Gl(0) — 
2{eT)'^ P2D /Tt-1 is obtained by after putting A = and fc = in Eq. (p7|). We see in Fig. 2 that the change from a 
Lorentzian and a Gaussian lineshapes depends also on the dimensionless ratio IAI/7. 

If we take the opposite limit to the inequality (p6|), i.e., we are dealing now with relatively short t^., then W{vpt/£c) 

in Eq. (|2^ ) has to be approximated by W{x)/de ~ yjTix~^ . As a result, the Lorentzian hneshape is obtained, from 



Eq. (PI), as 



G(A) = ^p,„-^ (29) 



2 

where the effective HWHM is now defined by 7^^.^ = 7[1 -I- (■\/7r/2)((5e ^c/S-wft)]- So, by increasing £c-, a transition 
from the Lorentzian to a Gaussian lineshape is obtained for be > 7. 

This peak modification is illustrated in Fig. 3, where G(A), calculated from Eq. (p5|), is displayed for fe/7 = 4.6. 
In Fig. 3, the solid, dashed, dotted, and dot-dashed curves correspond to hvp/^d = 15, 3.5, 0.7, and 0.2, respectively. 
Notice that the solid curve corresponds very closely to the Lorentzian given by Eq. (p9[), while the dotted and the 
dot-dashed curves are practically coincident Gaussians. In Fig. 4, we plot 7e///7, calculated from Eq. (|25|), as a 
function of hvp/id for decreasing values of 5e/^ = 4.6 (top), 2.3, 1.5, 1.1, 0.8, and 0.3 (bottom). We point out that 
all curves give 7e// in the local regime for Tivp/ic"! — 0. Thus, from Fig. 4, it is seen that nonlocal effects essentially 
make 7^// decrease for fe/7 ^ 1. 

Now we apply the present model calculation to interpret the experimental data of Ref. [||. We will consider the 
results in the low-temperature regime (less than 2 K), where the measured HWHM 7^^^ is practically independent 
of the temperature. We assume that only the inverted heterointerface for each QW has essential roughness due 
to one-monolayer variations (a « 2.5A) and, according to Appendix A, the characteristic energy is estimated as 
fe « 0.46 meV. The hard- wall model for a QW is used here to calculate e from ([Aq). For sample A, with electron 
density 1.6 x 10^^ cm~^, we assume that 7 ~ 0.1 meV (from mobility data), £c ~ 700 A and obtain from the pertinent 
curve, for fe/7 = 4.6, in Fig. 4, that 7e// ~ 0.22 meV, which coincides with the experimental data and corresponds to 
Tivf/£c1 = 15 (the solid square in Fig. 4). As a consequence, the Lorentzian lineshape given by the solid curve in Fig. 
3, is appropriate for sample A. Furthermore, for sample B, with density 1.5 x 10^^ cm~^, by assuming that 7 « 0.2 
meV and tc ~ 1400 A, we obtain 7,^^ « 0.45 meV from the pertinent curve, for fe/7 — 2.3, in Fig. 4, after using the 
calculated value hvpllc^ — 3.8 (the solid triangle). This value is in good agreement with the experimental result of 
Ref. For sample C, with density 0.8 x 10^^ cm~^ and assuming 7 k, 0.2 meV and Ic ~ 1000 A, we have 7^// ~ 0.45 
meV (indicated by the solid triangle in Fig. 4), i.e., the same as for sample B and also coincident with experimental 
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observations Then a good agreement with the experimental results is found for the case of single-side variations 
of heterointerfaces (see Ref. |16 about the case of two-side variations). 

Notice that the change from a local regime of tunneling (for long-range fluctuations of QW widths) to the general 
nonlocal case may be found by varying the temperature, that controls the relative contributions of homogeneous and 
inhomogeneous broadening^ In Fig. 5 we plot the lineshapes G(A), for temperatures 8 in the range 0.7 K - 10 
K, calculated from Eq. (psj), for sample B parameters taken from Ref. Now we have to add the thermal e - e 
scattering contribution to 7, leading to a renormalized value 7(0), which we approximate in the same way as in Ref. 

(see the solid curve in Fig. 3 of Ref. |^). In Fig. 5, the solid, dashed, dotted and dot-dashed curves correspond to 
8 = 0.7, 5, 7, and 10 K, respectively. One can see that by decreasing the temperature, the linewidth becomes smaller 
and the shape of peak is changed. While for 8 = 0.7 K, we observe a Lorentzian form of the peak due to strong 
nonlocal effects, manifested by inhomogeneous broadening induced by non-screened large-scale fluctuations, nonlocal 
effects are quite weak for 8 = 10 K. At this temperature the local regime prevails and the lineshape is the interplay 
of Gaussian and Lorentzian forms, given by Eq. (p7|), because in this case fc/7(8) w 0.5 and hvp/icliQ) ~ 0.8 the 
peak behavior is slightly more Lorentzian than Gaussian. 



IV. CONCLUDING REMARKS 



In the present work, wc have introduced a new electron scattering mechanism by long-scale non-screened roughness 
of heterointerfaces which contributes to the inhomogeneous broadening of the tunneling conductance peak in coupled 
DQWs. We have done a systematic analysis of this peak shape by taking into account the interplay between the 
introduced mechanism, and the usual homogeneous scattering broadening. A detailed comparison of the HWHM of 
the peak with experimental results revealed that the considered mechanism is relevant to interpret the data of Ref. 
1^, because in the experimental conditions strong nonlocal effects are manifested, through the proposed mechanism 
of inhomogeneous broadening, which modify drastically the lineshape (from the Gaussian to a Lorentzian) and the 
HWHM of the peaks. We call the attention to general considerations for the mechanism of appearance of non-screened 
long wavelength variations of the scattering potential, given in Appendix A, which should be relevant not only for 
the problem of inter-QW tunneling current, addressed here, but also for the study of general transport and optical 
properties of doped DQWs. 

Let us discuss the approximations used in our treatment. The single-electron approximation for the tunneling 
Hamiltonian in Sec. II is a generally accepted model and the expression for the tunneling current in the homogeneous 
case corresponds to the Bardeen's approach. ||l^ Here smooth variations of boundaries lead to changes of the levels 
of /— and r— QWs in the tunneling Hamiltonian, and we assume that small modifications of the tunneling matrix 
element can be neglected. For a discussion of the latter approximation, see Ref. The approximation for con- 
sidering the screening "on average" in the introduced large-scale potential consists in supposing that the correlation 
length ic is large in comparison with the Bohr radius and with transverse dimensions of the DQW structure as well. 
Since G(A) depends on the sum of the scattering broadening of different levels, we believe that the introduction of 
phenomenological parameters ^i^r instead of a detailed consideration the self-energy functions, does not lead to the 
omission of important contributions. We have also used the quasi-classical description for longitudinal motion which 
is valid when Sp ^ 7^// in the calculation of the path integrals in Eq. (|l^) and for integration over u and t in Eq. 
( ^ ) . We have assumed that variations of potential are sufficiently weak such that the acceleration of an electron due 
to long-scale random force (quasi-electric field) on a length of the order of £c is insignificant. We have considered 
that the inverted AlGaAs-GaAs heterointerface is much more rougher than the normal GaAs-AlGaAs heterointerface 
based on the results of Ref. jl^ for QWs similar to those used in DQW structures of Ref. 

To conclude, we now discuss the possibility of a more reliable test of the described mechanism of inhomogeneous 
broadening. In further experimental studies of the tunneling conductance it is necessary to make a more detailed 
analysis of the lineshape transition (from the Lorentzian to Gaussian form) , and comparison between HWHM data 
measured in samples grown in different conditions. Because the considered mechanism modifies also in-plane transport 
coefficients and optical properties of DQWs with long-wavelength inhomogeneities, further measurements as well 
theoretical studies of these phenomena are necessary. 

We believe that the present work establishes an essential contribution of large-scale non-screened fluctuations to 
the broadening of the tunneling conductance peak in DQWs in agreement with experimental results. 
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APPENDIX A: NON-SCREENED VARIATIONS OF THE RANDOM POTENTIAL 



Below we evaluate the non-screened random contributions to the potential ±5Ux that appears in Eqs. @ - (Q). 
The 2D Fourier transform of the screening potential, Vqz , is determined by the Poisson equation 



(- 



47re" 



(Al) 



where Srif^z is the concentrat ion induced by the total large-scale potential. Neglecting the overlap of I- and r- orbitals 
(pjz {j = I, r) we use in Eq. ( |Al[ ) the expansion 



j=l,r 



(A2) 



where Sn^j = —p2D(Uqj + V^j) is the in-plane induced concentration in the jth QW due to slow variations of the 
potential while the non-diagonal components of Sfix are small due to weak overlap of /- and r-orbitals. Note that, for 
the general case of slowly varying heterointerfaces, (pj^z should depend on the in-plane coordinate x but for the 2D 
case {sp <C Sj) such in -pla ne variations of orbitals are not important 
The solution of Eq.(Al) assumes the form 



Vqz - 



qz' I 



(A3) 



where for the case q ^ d ^ we have exp(— g|z — z'\) ~ 1. Then the diagonal components Vqj, j ~ l,r, are 

expressed in terms of the total concentration as 



dzip'jzVqz 



(A4) 



Substituting Suqj in Eq. ([A4|), we have 



V 



qj 



'—J2(Uc,s + Vq,), 

s—l,r 



(A5) 



where Or is the Bohr radius. Since the right-hand side of Eq. (A_5) does not depend on j we obtain V^i.r — Vq and 



we deal with an averaged screened potential Vq across the DQWs. If we assume large-scale variations, in particular 
(■c ^ Ois, and the condition qas/S <C 1, we derive the Eq. @), Vq ~ Uqj)/'2.. 

Now we present explicit expressions for the large-scale addends to the matrix Hamiltonian (|l]) due to the random 
variations of the DQWs heterointerfaces. Statistically independent boundaries variations are described by random 
functions 5^ and 5^ (see Fig. Ih) so that large-scale potentials take the form C/j^rx — 2ei^r(^x ~ ^^)/di,r, where 
ei,r are the energies of levels in Z, and r-QWs with width di ^- As it was shown in Ref. 1^, the weak contributions 
due to variations of the tunneling matrix element may be neglected in the evaluation of Eq. (^ up to second-order 
in T . Substituting C//..rx into the correlation functions in Eq. (^) and supposing that all interfaces are statistically 
equivalents, we obtain the correlation functions as 



2ea 



exp 



(A6) 
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where e ~ r, ■ 



is the correlation length. Finally, 

2 /«2l 



ii,r, and a is the averaged deviation of heterointerfaces and 
the correlation function w(|...|) in Eq. (Q) takes the form w(|x — x'|) — (fc) exp[— (x — x')^ /i'^] with a characteristic 
energy Se = 4e a/d for the case of two-side variations. However, 5e = 2^/2e a/d for the o ne-si de variation case in a 
QW when or is equal to zero and one of the two correlation functions given by Eq. (A6) is set to zero. 



APPENDIX B: OPTIMAL FLUCTUATION METHOD 



The evaluation of the Eq. (2^) is based on the separation of the optimal trajectory (with maximal contribution 
to the path integral) and on the comparison of typical variations to such a trajectory, (5x^ and Sy-^ with £c. First 
variations of the actions in the exponential factors of Eqs. (|l8|), (|l^) with respect to path variations Sx^ and (5y^ are 
written as follows 



2ne 







r dr' 


H 


'o 





2Ml 



dr' [u(t - T')/ti + x^ - x^'] w{\u{t - r')/ti -|- x^- - x^'|) 



u 



(1±1) 



wi\uT/ti + X, - ur'/ia - y^' - - (1 ± 1) |)}, 



and 



(Bl) 



drSyr ■ {±myr - 



2Ul J, 



dT'[±u(T - T')/t2 +yr- yAw{\ ± u(t - T')/t2 + Yr - Yr' 



2Ul j„ 



dr' 



t2 h 



w{\nT'/h + ^r' - UT/t2 - + - (1 ± 1) I)}, 



(B2) 



where the upper (lower) signs correspond to upper (lower) signs in Eqs. (jig), (^^ and we have transformed correlation 
functions as follows 



w(|x + 5x|) - u)(|x|) ~ -2x • 5x w{x)/el. 



(B3) 



Since (Jx^. and Sy^ are arbitrary variations, the optimal trajectories are determined by the system of Euler-Lagrange 
equations as 



2MI J, 



C?T'[u(r - T')/ti +Xr - X^']w(|u(t - t')/<i + X^ - Xr> |) 



and 
± my^ 



/ ' rfT'[u(T/ti - r'/ia) + X, - y,, - u (1 ± 1) /2]w{\uT/h + - ur'/ts - y^' - u (1 ± 1) /2|), (B4) 



dT'[±u(T - r')/t2 + yr - yr']w{\ ± u(t - T')/t2 + Y. " Yr'l) 
tl 

dT'[u{T/t2 ~ r'/h) -xr'+yr-u{l± 1) /2]w{\uT/t2 + Yr - ur'/ii - X,, - u (1 ± 1) /2|). (B5) 



i 



2MI Jo 



In order to estimate the maximal deviations, imax, Umax, we suppose that UT^2^/ti^2 ^ |a^max|, |2/max| for typical 
parameters used in calculating G(A) such that Xmax = a^rj""" and j/max = J/tJ""=- Thus, the right-hand sides of Eqs. 
(B4), (B5) do not depend on x^ and y^ and they can be easily integrated with the boundary conditions x^-^o,*! = 



-^2 



and Yr=o.t2 = 0. For the upper-sign contributions, estimating uw(|u|)| « £cS£ we transform the right-hand sides of 



Eqs. (B4), (B5) into iSe {\ti\ + \t2\)/{2mh£c)- Thus, the result for the maximal deviations is 



'^max 
ymax 



Se {\ti\ + \t2\) 



(B6) 



This satisfies the condition |a;max|, |ymax| ^ £c for the upper-sign contributions because ti^2 are estimated as h/Sf 
Indeed, then we have |a;max|/^c and |ymax|/^c ~ {Se/ep)'^ x (£c/4ei?) <C 1, where the characteristic energy Sc ~ 
{h/£cf/2m < e^. 
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A more c aref ul consideration is necessary for the lower-sign contributions when \ti — ^ Ji/sp, such that ti.2 
Thus, Eq. (|B4| ) can be rewritten as 



2 

iuSe t 



dz'{z - z') exp [-u^{z' - zf/e'; 



(B7) 



where z = r/t. Here, from Eq. (B5), it follows that ~ 0, and taking into account the boundary conditions, 
yT=o,t2 = 0; we are led to yr ~ 0. It is easy to solve Eq. (B7) taking into account in its right hand side that for the 
lower-sign contributions i r, as it was shown in Sec. IV, due to the fact that ^ 7^//- Then solving Eq. (|B7|), 
with boundary conditions x^=o,ti = Oi we obtain the maximal deviation |xmax| as 



IGmhi. 



(B8) 



where Fx{z) = [1 — exp(— z^)]/z and Fx{z) « z, for <C 1, while Fx{z) « 1/z, for z^ ^ 1 and the maximum of 
F^{z) < 0.65 corresponds to z close to the unity. This estimation practically coincides with Eq. (B6) if u/£c ~ 1, 
while the maximal deviations are smaller then the results in Eq. (B6) both for u/£c ^ 1 and u/ic ^ 1- 



°0n leave from: Institute of Semiconductor Physics, Kiev, National Academy of Sciences of Ukraine, 252650, 
Ukraine 



R. Tsu and L. Esaki, Appl. Phys. Lett. 22, 562 (1973). 

P. Mazumder, A. Kulkarni, M. Bhattacharya, J. P. Sun, and G. I. Haddad, Proc. IEEE 86, 664 (1998). 
J. Smoliner, W. Demmerle, G. Berthold, E. Gornik, G. Weimann, and W. Schlapp, Phys. Rev. Lett. 63, 2116 (1989). 
J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Appl. Phys. Lett. 57, 2324 (1990). 
S. Q. Murphy, J. P.Eisenstein, L. N. Pfeiffer and K. W. West, Phys. Rev. B 52, 14825 (1995). 

J. A. Simmons, M. A. Blount, J. S. Moon, S. K. Lyo, W. E. Baca, J. R. Wendt, J. L. Reno, and M. J. Hafich, J. Appl. 
Phys. 84, 5626 (1998). 

F. T. Vasko, J. P. Sun, G. I. Haddad, and V. V. Mitin, J. AppL Phys. 87, 3582 (2000). 
O. G. Balev, F. T. Vasko, F. Aristone, and N. Studart, submitted to Phys. Rev. B. 

F. T. Vasko and O.E. Raichev, Phys. Rev. B 50, 12195 (1994); F. T. Vasko and O. E. Raichev, JETP 80, 539 (1995). 
L. Zheng and A. H. MacDonald, Phys. Rev. B 47, 10619 (1993). 

O. E. Raichev and F. T. Vasko, J. Phys.: Cond. Matter, 8, 1041 (1996); ibid. 9, 1547 (1997). 
L. S. Schulman, Techniques and Applications of Path Integration, John Wiley & Sons, Inc., New York, (1981). 
B. I. Halperin and M. Lax, Phys. Rev. 148, 722 (1966). 

Notice that the general form of Eq. (|2^) may be obtained by averaging the Lorentz peak with a random level splitting Ax. 
The essential point here is the condition of the local regime, given by Eq. (p^). 
L. N. Pfeiffer, E. F. Schubert, and K. W. West, Appl. Phys. Lett. 58, 2258 (1991). 

For two-side variations of heterointerfaces, as in Eq. (AC), we obtain £c = 350 A, 700 A and 500 A for samples A, B and 
C respectively. Because in Ref. it is reported that 200 A QWs are separated by the tunnel barriers with 175 A to 350 
A (even though the real values are not given samples for A, B and C) and considering the results of Ref. [^, we have 
assumed here only one-side fluctuations of the QW for samples in Ref. 
[17] J. Bardeen, Phys. Rev. Lett. 6, 57 (1961). 



10 



FIGURE CAPTIONS 



Fig. 1. a) Spatial variations of the energy levels in I- and r-QWs along the x direction without screening (dotted 
curves) and with screening (solid curves); b) Band diagram of DQWs, along the z direction, with nonideal heteroint- 
erfaces shown by dashed lines. 

Fig. 2. Modified lineshapes G(A), taken from Eq. (p7|), normalized by Gl(0) = 2(eT)^p2_D/?i7, when nonlocal 
effects are negligible, for different contributions of short-range scattering, characterized by the phenomenological 
broadening parameter 7, and non-screened large-scale disorder, characterized by Ss and ic- The solid, dashed, dotted, 
dot-dashed, and dot-dot-dashed curves correspond to Se/^ = 0.3, 0.6, 1, 3, and 6, respectively. The solid curve is 
almost a Lorentzian lineshape, while the dot-dot-dashed curve is very close a Gaussian curve. 



Fig. 3. Modified lineshape G(A), calculated from Eq. (pSD, when nonlocal effects are taken into account, for 
fc/7 = 4.6 and different values of hvp/tc'^- The solid, dashed, dotted, and dot-dashed curves correspond to hvp/ic'^ = 
15, 3.5, 0.7, and 0.2, respectively and represents, for instance, the increase of £c- Note that the lineshape evolves 
practically from a Lorentzian (solid curve), given by Eq. (^), to a Gaussian one (dot-dashed curve). 



Fig. 4. Dimensionless tunnel resonance width 7e///7 as function of hvp/f^cj, calculated from Eq. (|25|). The solid 
curves from top to bottom correspond to Se/j = 4.6, 2.3, 1.5, 1.1, 0.8, and 0.3. For hvp/^cl = 0, nonlocal effects are 
negligible. 



Fig. 5. Modified lineshape G(A), for data of sample B of Ref. calculated from Eq. (25), in which temperature 
effects were included. The solid, dashed, dotted and dot-dashed curves correspond to temperatures 0.7, 5, 7, and 10 
K, respectively. In G_l(0) it was used 7 = 0.2 meV corresponding to the sohd curve. 
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